Fluctuations and Rheology in Active Bacterial Suspensions 
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We probe non-equilibrium properties of an active bacterial bath through measurements of corre- 
lations of passive tracer particles and the response function of a driven, optically trapped tracer. 
These measurements demonstrate violation of the fluctuation-dissipation theorem and enable us to 
extract the power spectrum of the active stress fluctuations. In some cases, we observe scaling 
in the noise spectrum which we show can be derived from a theoretical model incorporating coupled 
stress, orientation, and concentration fluctuations of the bacteria. 

PACS numbers: 82.70.Dd, 83.60.Bc, 5.40.-a,87.17. Jj,87.18.Hf 



Active complex fluid systems such as living cells [ll, [3], 
assemblies of motors and filaments flocks of birds 
and vibrated granular media Q differ from conven- 
tional equilibrium media in that some of their compo- 
nents consume and dissipate energy, thereby creating a 
state that is far from equilibrium. An understanding of 
model active systems, even at a phenomenological level, 
provides insight about fundamental non-equilibrium sta- 
tistical physics and, potentially, about the inner workings 
of biological systems. Bacterial baths [1, 0, H, H, US] are 
attractive model active systems because energy input is 
homogeneous, because individual bacteria can be directly 
observed, and because critical parameters such as density, 
activity, and swimming behavior Q can be brought under 
experimental control. Indeed, recent experiments have 
reported on a rich variety of non-equilibrium phenom- 
ena in this system class including anomalous diffusion 
01 and pattern formation [llj, while recent theories 
of self-propelled organisms predict ordered phases such 
as "flocks" their instabilities [l3|, novel rheological 
effects [l^, and giant density fluctuations p^ . 

In this Letter, we describe measurements of the fluc- 
tuations and mechanical responses of an active bacterial 
suspension. In contrast to previous work 0, H, 0, we 
concurrently measure the one- and two-point correlation 
functions of embedded passive tracer particles to assess 
material fluctuations over a wide range of length scales 
p^ . and we independently measure the macroscopic re- 
sponse function of the "active" medium by displacing a 
particle through it with an optical trap JGJ. Even at a 
low volume fraction ((/) ^ 10"'^) of bacteria, fluctuations 
in the medium are substantially greater than they are in 
the absence of bacteria while rheological response is un- 
changed, implying a strong violation of the fluctuation- 
dissipation theorem (FDT). The mean-square displace- 
ments (MSDs) of tracer particles as a function of time 
At, depend on swimming behavior. For wild-type bacte- 
ria, the MSD extracted from two-point correlations scale 
superdiffusively as At"^^^ for the time scale of our experi- 
ments, and the stress power spectrum 2] as a function of 
frequency uj scales as (jj/^ytJ. Existing theories of active 
media [121 predict long-time tails and anomalous correc- 
tions to diffusion but not superdiffusion. We propose a 
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FIG. 1: 1-pt (open symbols) and 2-pt (closed symbols) mean 
square displacements divided by the lag time. At, for 2a 
= 2 fim particles in bath of wild-type RP437 (triangles) 
and tumbling RP1616 (circles) bacteria at (p — .003. In- 
set: RDrr{R, At = .067 sec) demonstrating Drr ~ 1/R for 
i? > 10 nm, and implying correlation length < 10 /im. 



theoretical model, following Ref. [13], that accounts for 
our experimental observations. 

We used two strains of E. coli, a rod-shaped bac- 
terium with dimensions 3x1 /im, in these studies: 
RP437, the "wild-type", which runs and tumbles (l7j 
and RP1616, the "tumbler", which predominantly tum- 
bles [31 . Overnight cultures were diluted 1/300 in Luria 
Broth (Difco) and grown at 25 °C for 6 hrs. Subsequently, 
they were centrifuged for 10 minutes at 5000 rpm and re- 
suspended to the desired concentration in a buffer com- 
prised of 10 mM K2HPO4, 0.1 mM EDTA, and 0.2 wt 
% glucose (pH = 8.2), which was added to maintain bac- 
terial motility. We added a small amount {(f>s = 10~'^) 
of fluorescently labeled polystyrene spheres (Duke Sci- 
entific) of radius a to the bacterial suspension, and to 
density match them, we added 15 wt % sucrose to the 
solution. To prevent bacterial adhesion, we prepared the 
chambers from BSA coated glass slides and coverslips. 
We used parafilm spacers to bring the thickness of the 
chambers to ~ 240 /im; and we recorded quasi-2D image 
slices from the middle of the 3D chamber. Samples were 
loaded into the chamber and sealed with optical glue just 
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prior to each run. 

We quantified the fluctuations in the bacterial bath by 
computing MSDs from the motions of embedded micron- 
sized tracers [15|]. The one-point displacement (MSDl) 
is defined by {Ar^At))^ = ^^(Ar,(t, At)Ar,(t, Ai)), 
where Ari{t,At) = ri{t + At) — ri(t) is the par- 
ticle displacement in the i = {x, y, z) direction in 
time Ai, and the brackets denote time and ensem- 
ble averaging. The two-point displacement (MSD2) 
is defined as (Ar2(At))2 = (2R/a) Drr{R, At), where 
Drr{R, At) is the longitudinal component Drr{R, At) — 
DijRiRj/R^ of the two-point tensor Dij{R,At) = 

( Arf ^ (t, At)Arf ^ {t, At)), which measures correlations 
of two distinct particles (1,2) with an initial separation 
R. Over the time scale of our experiments, R lies in the 
focus plane of our microscope and its magnitude R = |R| 
is greater than that of individual particles' displacements. 
The main advantage of two-point microrheology is that 
it provides a more reliable measure of length scale de- 
pendent fluctuations in media where the length scale of 
heterogeneities and tracer boundary conditions are not a 
priori known 0, [l5| . Indeed, since Dij {R, At) is ensem- 
ble averaged over tracer pairs with R^ a, it reflects the 
dynamics of the medium on larger length scales than the 
tracer size, permitting quantitative measurements even 
in the presence of heterogeneities. In general, MSD2 will 
equal MSDl if heterogeneities in the medium have length 
scales smaller than the tracer size, otherwise they will 
differ in both magnitude and functional form. 

Typical MSD data are presented in Fig.[l] which shows 
that the one-point MSD in both bacterial strains dis- 
plays a crossover from superdiffusive behavior at short 
lag times ((Ar^)i ~ Ai",l < a < 1.5) to diffusive be- 
havior (a = 1) at long lag times. This observation is 
similar to that of Ref. [7]. However, our two-point data 
for the wild-type, by contrast, exhibit a nearly power-law 
superdiffusion ((Ar^)2 ~ At^-^) over 2.5 decades of ob- 
servation time. We also verified that Drr{R, At) ^ 1/R 
(see inset of Fig. [1]) , indicating that the bacterial bath, 
though an active medium, can be treated on the separa- 
tion scale i? as a coarse-grained continuum whose proper- 
ties can be probed with two-point microrheology 0, [Tst . 
That MSDs exhibit superdiffusion is suggestive of but not 
a proof of violation of FDT, which requires an indepen- 
dent measurement of the rheological response function. 

Response measurements were performed using an os- 
cillating optical tweezer setup similar to that of Ref. 
[T^ . Briefly, an optical trap with typical trap stiffness 
of 1 X 10~^pN/nm was formed by focusing an ~ 100 
mW laser beam (A = 1054 nm) through a 1.3 NA oil im- 
mersion objective (Zeiss). The trapping beam position 
was sinusoidally scanned using a galvo-mirror at frequen- 
cies from 0.5 to 500 Hz. A 4.0 /im PS sphere was trapped 
~ 6 fim from the coverslip. The position of the tracer was 
detected using forward scattered light from a co-linearly 
aligned HeNe laser beam focused onto a split photodi- 
ode (Hamamatsu S4204). The photodiode signal was fed 
into a lock-in amplifier (Stanford Research Systems 530) 




0.1 1.0 10.0 100.0 1000.0 

to (rad/s) 



FIG. 2: (color online) Frequency dependent viscosity derived 
from oscillating trap measurements for 4.0 nm sphere in wa- 
ter (solid squares), the tumbler (solid circles), and the wild- 
type (solid triangles) at = .003. Viscosities r]2{u!) derived 
from the averaged two-point measurements using the general- 
ized Stokes-Einstein relation are plotted for the tumbler (open 
circles), the wild- type (open triangles), and a bead in water 
(open squares). Inset: Normalized Displacement of a 4.0 fim 
sphere in the optical trap as a function of driving frequency for 
wild-type (triangles), tumbler (circles), and water (squares). 
Line is a fit to d(tj) (see text). 



along with the reference from the driving function gener- 
ator signal. The displacement and phase of the trapped 
particle output by the lock-in amplifier were logged into 
a PC running Lab View (National Instruments). 

The equation of motion for a particle of radius a 
trapped in an oscillating harmonic potential may be writ- 
ten as: Gnrjax = ~k [x — Acos{u!t)], where ry is the vis- 
cosity of the medium, k is the stiffness of the trap, and A 
is the driving amplitude. Its steady state solution yields 
the normalized displacement of the sphere in the trap: 

d{ijj) = 1 + [6TTari{LLj)uj/k] > 

The inset of Fig. [5] shows the raw normalized dis- 
placement data for a particle in water and for a par- 
ticle in an active bacterial bath (RP1616). The solid 
line is a fit to d{uj) with rj = 0.001 Pa • s, trap stiffness 
k = 8x 10~^pN/nm, and radius a = 2.0 /xm. Both sets 
of experimental data agree with each other and with the 
theoretical curve. From them, we extract the viscosity 
ri{uj) shown in the main graph of Fig. [21 Clearly, the 
presence of actively swimming bacteria at volume frac- 
tion lO"'^ does not modify the viscosity of the medium 
from that of water, ri{uj) = rjo = 0.001 Pa • s. We mea- 
sured the same ?7(a;) using a capillary viscometer. 

While recent theories of active systems predict novel 
enhancement in the viscosity [isj . our experiments are 
well below the concentration at which these effects are 
observable. Instead, our results are consistent with 
the Einstein result for hard spheres: rj ~ rio{l + ^(j)), 
namely, a negligible modification in the viscosity for 
(f> 10"'^. Moreover, assuming for the moment the gener- 
alized Stokes-Einstein relation, we can extract the (FDT 
consistent) response from the collapsed two-point dis- 
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FIG. 3: (color online) (a) Collapsed 1-pt (open symbols) and 
radius collapsed averaged 2-pt (closed symbols) MSDs for the 
tumblers at <^ = .003. The solid line is the master curve: 
7+(l— 7)(1 — e^^ )/x. Circles, triangles, and squares are 
for particle diameters 2a = 2.0, 5.0, and 10.0 fim, respectively. 
Inset; Raw 1-pt (open symbols) and 2-pt (closed symbols) 
MSDs for tumblers, (b) Radius rescaled 2-pt (closed symbols) 
MSDs for the wild-types ai (f) = .003. Circles, triangles, and 
squares are for particle diameters 2a = 2.0, 5.0, and 10.0 fj.m, 
respectively. Inset: Radius rescaled 1-pt MSD (open symbols) 
for the same (f) and particle size data. 



placement (MSD2) [15]: r]2{uj) = kBT/3Truj'^a{Ar'^{uj))2, 
as shown in Fig. O The difference between r/(aj) and 
r?2(w) explicitly indicates a strong violation of FDT. 

Next, to access the heterogeneity of the bacterial bath, 
we explored the length-scale dependence of fluctuations 
by systematically varying the size of the tracers at a fixed 
bacterial concentration. The inset of Fig. [3^ shows MSDs 
obtained for spheres in the tumbler bath. All samples 
and all tracer sizes exhibit a crossover from superdif- 
fusion to diffusion on similar timescales, with an en- 
hanced diffusion coefficient D = 7I?t, where 7 — 4.3 and 
Dt — ksT / {QirriQa) is the equilibrium coefficient. More- 
over, MSDl and MSD2 are nearly equal in magnitude and 
functional form. Rescaling time by the crossover time 
r and the MSDs by 20x^1 collapses all the data onto 
a master curve: [{/^rl(M))i,{/\r'^{M))2\/{2DT/^t) = 
7 + (l — 7)(1 — e""^ )/x, where x — At/r. Figure [3^ 
shows the collapsed MSD data along with the master 
curve with r = 0.1 s. 

The MSDs for the wild- type are strikingly different: 



the MSDl exhibits a crossover dependent on tracer size, 
while all of the MSD2 exhibit superdiffusion with nearly 
the same exponent of 1.5 over 2.5 decades of time, 
independent of the tracer size. The trivial rescaling 
a{Ar^{At))2 coUapsed the respective MSD2 data [Fig. 
12)3]. Under this rescaling, however, (and other simple 
scaling functional forms as well) the wild- type MSDl 
failed to collapse [inset of Fig. [HJs], signaling the pres- 
ence of heterogeneity on the tracer length scale. The 
superdiffusive exponent of the MSDl approaches that of 
the two-point data (a ~ 1.5) as a increases. This sug- 
gests that one-point measurements are intrinsically am- 
biguous: the activity inferred depends on the tracer size 
and boundary conditions 0, [l^. Two-point measure- 
ments, in contrast, provide a more robust characteriza- 
tion of the long-wavelength fluctuations of the medium 
than one-point measurements. 

We employ a recently developed phcnomcnological the- 
oretical framework for an active medium to interpret the 
experimental MSD data . The bacterial activity gives 
rise to non-thermal stress fluctuations whose power spec- 
trum A(u;) can unambiguously be extracted from two- 
point microrheology data via 
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The results are exhibited in Fig. |4^. For water, the 
power spectrum is flat. For the tumblers, it is nearly 
Lorentzian, flat at low frequencies with a knee at higher 
frequencies. For wild-types, it exhibits power-law behav- 
ior, A(a;) ~ w""'^, over 2.5 decades. In both cases, A(a;) 
is substantially greater than in a thermal system. For the 
wild-type, the prefactor Ao of A(a;) rises linearly with the 
bacterial concentration, as shown in Fig. |4|d. 

We propose a simple model to account for the observed 
spectra of the tumblers and the wild-types within the the- 
oretical framework of Ref. [l^. In the process of tum- 
bling or swimming, each bacterium contributes an ad- 
ditional active stress to the medium. It has the form: 
a{j oc c{x,t)Sij{x,t), where c(x, t) is the concentration 
of the bacteria and Sij is a force-dipole density gener- 
ated by the active bacteria. In wild types and tumblers, 
forces are directed, respectively, along and perpendicu- 
lar to the long-bacterial axes. Thus, in wild-types, Sij 
is equal to the uniaxial nematic order parameter 
whereas in tumblers, it is equal to a biaxial order param- 
eter Q^. Active processes enhance stress fluctuations, 
and assuming long-range isotropy, the active stress fluc- 
tuations can be expressed as ( cr^(q, a;)cr^;(— q, — a;) ) = 
A^(g, Lo) [6ik6ji + SiiSjk - | for both tumblers and 

wild- types. The power spectrum in Eq. ([1]) is related to 
A„{q, cj) by A{uj) = A<,(q = 0, uj). 

The linearized equation for Qf, {A = U, B) is dtQ^ — 



where ta is the relaxation 



time, the correlation length of Q^-, and Sij is a spatial- 
temporal white noise with zero mean. Interactions 
among bacteria favor long-range order in QY^ but not in 
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FIG. 4: (a) The spectrum A(a;) of active stress fluctuations 
obtained from two-point microrheology and active response 
measurements. The triangles are the wild-types, circles are 
the tumblers (both — .003), squares are water (0 = 0). (b) 
Linear dependence of the prefactor Ao in A{ui) on the volume 
fraction (f> of the wild-type bacteria; At = 2irjkBT. 



Qfj, implying that tb <C tjj and £,b ^ £,u- In both cases, 
the concentration of bacteria obeys the continuity equa- 
tion: dtSc = — V- J with Ji = —DdiSc — a2CodjQfj +SJi, 
where cq is the average concentration, D is the diffu- 
sion constant, SJi is a random current, and the second 
term stems from the nonequilibrium driving of mass flow 
These equations lead after Gaussian decoupling to 



Ao-(q,w) = Acr{q^crA,^^) with ^^tA ^ Ca- In tumblers 
is very small, and Aa-iu) can be replaced by a Lorentzian 
A(aj) with characteristic time tb ^ 0.1s (Fig. |4]) [23| 
in both MSDl and MSD2 implying agreement between 
{AP{At))i and {AP{At))2. In wild- types, A^iq^u,^) 
can be replaced by A{u!) in MSD2 when R > [with 
< 10 /im (see inset of Fig. [1])], but not in MSDl when 
it probes lengths shorter than ^t/, implying different val- 
ues for (Ar2(At))i and (Ar^{At))2 [Fig.[l]. Our calcula- 
tions yield A(uj) ^ cq/^/lu for ujtij > 1 to lowest order in 
Co, in agreement with our experiments (20| . This result 
arises essentially from the concentration fluctuations in 
the active stress, which were ignored in previous theories. 
Note also that one might expect that swimming bacteria 
have a tendency to develop long-range polar rather than 
the nematic order of our model, but we find that incipient 



polar order yields A(lu) ^ Cg < 
with our measurements [20l] . 
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in clear disagreement 



In conclusion, using a combination of passive 2-pt mi- 
crorheology and active response measurements, we have 
demonstrated that the macroscopic stress fluctuations 
depend sensitively on microscopic swimming behavior of 
the bacteria. When contrasted with other active systems 
such as living cells, our results suggest that departures 
from equilibrium proceed via non-universal mechanisms. 
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